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Abstract 



A sample of 3120 K± — ¥ TT ± fi + /i decay candidates with (3.3 ± 0.7)% background 
contamination has been collected by the NA48/2 experiment at the CERN SPS, allow- 
ing a detailed study of the decay properties. The branching ratio was measured to be 
BR = (9.62 ± 0.25) x 10~ 8 . The form factor W(z), where z = {M^/M K f, was pa- 
rameterized according to several models. In particular, the slope of the linear form factor 
W{z) = W {1 + Sz) was measured to be S = 3.11 ± 0.57. Upper limits of 2.9 x 10~ 2 and 
2.3 x 10 -2 on possible charge asymmetry and forward-backward asymmetry were established 
at 90% CL. An upper limit BR(K ± -> ir^^^) < 1.1 x 10" 9 was established at 90% CL 
for the rate of the lepton number violating decay. 
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Introduction 

The flavour-changing neutral current decays — > n £ + £~ (denoted K^m below, £ = e, /i), 
induced at the one-loop level in the Standard Model (SM), are well suited to explore its structure 
and, possibly, its extensions. The rates of these transitions are dominated by the long-distance 
contributions involving one photon exchange. They have been described in the framework of 
Chiral Perturbation Theory (ChPT) [1] in terms of a vector interaction form factor (which 
characterizes the dilepton invariant mass spectrum) determined by experimental measurements. 
Several models for form factor have been proposed 0[3j|4]. 

The first experimental observation of the K£„„ process was published by the BNL E787 
collaboration in 1997 |5j. It was followed by a BNL E865 measurement |6j which established 
the vector nature of the decay, and found the form factor and the decay rate to be in agreement 
with the expectation based on the earlier K nee measurements 018]. The most stringent upper 
limit on the rate of the lepton number violating K + — > tt~ fx + fi + decay also comes from the 
E865 experiment [9]. Later the HyperCP experiment analysed samples of both and K~„ 

decays |10| . which, in addition to decay rate measurements, allowed setting a limit on the CP 
violating rate asymmetry. The total K^un sample collected by the three experiments amounts 
to ~ 700 candidates. 

A new measurement of the K^„„ decay based on the data collected by the NA48/2 experiment 
at the CERN SPS in 2003-2004 is reported in this letter. The event sample is ~ 4.5 times larger 
than the total world sample, and has low background contamination, allowing form factor, rate 
and asymmetry measurements at an improved precision. 

1 The NA48/2 experiment 

The NA48/2 experiment, specifically designed for charge asymmetry measurements [llj, uses 
simultaneous K + and K~ beams produced by 400 GeV/c primary SPS protons impinging on a 
beryllium target. Beam particles with momentum (60 ± 3) GeV/c (r.m.s.) are selected by an 
achromatic system of four dipole magnets with zero total deflection ('achromat'), which splits 
the two beams in the vertical plane and then recombines them on a common axis. The beams 
pass through momentum defining collimators and a series of four quadrupoles designed to focus 
the beams at the detector entrance plane. Finally the two beams are again split in the vertical 
plane and recombined in a second achromat. 

The beams enter the fiducial decay volume housed in a 114 m long cylindrical vacuum tank 
with a diameter of 1.92 m upstream, increasing to 2.4 m downstream. Both beams follow the 
same path in the decay volume: their axes coincide within 1 mm, while the transverse size of 
the beams is about 1 cm. With 7 x 10 11 protons incident on the target per SPS spill of about 
4.8 s duration, the positive (negative) beam flux at the entrance of the decay volume is 3.8 x 10 7 
(2.6 x 10 7 ) particles per pulse, of which 5.7% (4.9%) are K + (K~). The K+/R- flux ratio is 
1.79. The fraction of beam kaons decaying in the vacuum tank at nominal momentum is 22%. 

A detailed description of the NA48 detector can be found in Ref. [12]. The decay volume is 
followed by a magnetic spectrometer housed in a tank filled with helium at nearly atmospheric 
pressure, separated from the vacuum tank by a thin (~ 0.4%Ao) Kevlar® window. A thin- 
walled aluminium beam pipe of 16 cm outer diameter traversing the centre of the spectrometer 
(and all the following detector elements) allows the undecayed beam particles and the muon 
halo from decays of beam pions to continue their path in vacuum. The spectrometer consists 
of four octagonal drift chambers (DCH) composed of eight planes of sense wires: DCH1, DCH2 
located upstream, and DCH3, DCH4 downstream of a dipole magnet. The magnet provides 
a horizontal transverse momentum kick Ap = 120 MeV/c for charged particles. The spatial 
resolution of each DCH is a x = a y = 90 /im. The momentum resolution of the spectrometer is 
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a p / P = (1.02 0.044 • p)% (p in GeV/c). 

A plastic scintillator hodoscope (HOD) used to produce fast trigger signals and to provide 
precise time measurements of charged particles is placed after the spectrometer. The HOD has 
a regular octagonal shape, and consists of a plane of vertical strip-shaped counters followed by 
a plane of horizontal ones (128 counters in total). 

The HOD is followed by a liquid krypton electromagnetic calorimeter (LKr) used for particle 
identification in the present analysis. It is an almost homogeneous ionization chamber with an 
active volume of 7 m 3 of liquid krypton, segmented transversally into 13248 projective cells, 
approximately 2x2 cm 2 each, 27 Xq deep and with no longitudinal segmentation. The transverse 
sizes of the HOD, DCHs and LKr are about 2.4 m. 

A muon detector (MUV) essential for muon identification in the present analysis is located 
further downstream. The MUV is composed of three planes of plastic scintillator strips (aligned 
horizontally in the first and last planes, and vertically in the middle plane) read out by pho- 
tomultipliers at both ends. Each strip is 2.7 m long and 1 cm thick. The widths of the strips 
are 25 cm in the first two planes, and 45 cm in the third plane. The MUV is also preceded 
by a hadronic calorimeter (not used for the present measurement), which is an iron-scintillator 
sandwich with a total iron thickness of 1.2 m. Each MUV plane is preceded by an additional 
0.8 m thick iron absorber. 

A dedicated two- level trigger has been designed for collection of three-track decays. A 
description of the trigger algorithm and the sources of its inefficiency, which is typically at the 
10~ 3 level and is therefore negligible for the present analysis, can be found in Ref. |llj . 

A detailed GEANT3-based [13] Monte Carlo (MC) simulation which includes full detector 
geometry and material description, stray magnetic fields, DCH local inefficiencies and misalign- 
ment, detailed simulation of the kaon beam line, and time variations of the above throughout 
the running period is used to compute the acceptances for signal, normalisation, and background 
channels. 

2 Event selection 

The K n nu rate is measured relative to the abundant — > ■K ± it Jr ir~ normalisation channel 
(denoted below). The K nfJifl and samples are collected concurrently using the same 
trigger logic. The fact that the ^ and tt^ masses are close {m^/m^ = 0.76) results in similar 
topologies of the signal and normalisation final states. This leads to first order cancellation of 
the systematic effects induced by imperfect kaon beam description, local detector inefficiencies, 
and trigger inefficiency. 

Selection conditions 

Three-track vertices (compatible with either K nfJifl or decay topology) are reconstructed 
by extrapolation of track segments from the spectrometer upstream into the decay volume, 
taking into account the measured Earth's magnetic field, stray fields due to magnetization of 
the vacuum tank, and multiple scattering. 

The selection procedures for the K nfJjfJi and K^ n modes have a large common part: the 
presence of a vertex satisfying the following criteria is required. 

• The vertex longitudinal position is within the fiducial decay volume (i.e. downstream the 
final collimator). 

• The vertex tracks are required to be consistent in time (within a 10 ns time window), 
consistent with the trigger time, and to be in DCH, HOD, LKr and MUV geometric 
acceptances. Track momenta are required to be above 10 GeV/c to ensure high muon 
identification efficiency. Track separations are required to exceed 2 cm in the DCH1 plane 
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to suppress photon conversions, and 20 cm in the LKr front plane to minimize particle 
misidentification due to shower overlaps. 

• The total charge of the three tracks is Q = ±1. 

• The total momentum of the three tracks | Y^,Pi\ 1S consistent with the beam nominal range: 
(54; 66) GeV/c. 

• The total transverse momentum of the three tracks with respect to the mean beam direc- 
tion (which is precisely measured using the sample) is p\ < 0.5 x 10~ 3 (GeV/c) 2 . 

If several vertices satisfy the above conditions, the one with the lowest fit x 2 is considered. 
The candidates are then selected using the following particle identification and kinematic 

criteria. 

• The vertex is required to be composed of one ir^ candidate (with the ratio of energy 
deposition in the LKr calorimeter to momentum measured by the spectrometer E/p < 0.95, 
which suppresses electrons, and no in-time associated hits in the MUV), and a pair of 
oppositely charged ^ candidates (with E/p < 0.2 and associated hits in the first two 
planes of the MUV). The muon identification efficiency has been measured to be above 
98% for p > 10 GeV/c, and above 99% for p > 15 GeV/c. 

• The invariant mass of the three tracks in the 7r ± ^ + /i~ hypothesis lies in the range {M^^ — 
Mk\ < 8 MeV/c 2 , where Mx is the nominal charged kaon mass. 

Independently, the following criteria are applied to select the K^ n sample. 

• The pion identification criterion described above is applied to a single pion only, to sym- 
metrize the selection of the signal and normalisation modes and diminish the corresponding 
systematic uncertainties. 

• The invariant mass of the three tracks in the 3^ hypothesis lies in the range \M^—Mk\ < 
8 MeV/c 2 . 

No restrictions are applied to the additional energy deposition in the LKr calorimeter, which 
decreases sensitivity to accidental activity. 

Signal sample 

The reconstructed 7r ± /x + /i~ invariant mass spectrum is presented in Fig. [1^,: a K^aa decay signal 
is observed. The number of K ntlfl candidates in the signal region is = 3120, of which 2003 

(1117) are K + (K~) candidates. The measured M wflll resolution is a nllfl = 2.5 MeV/c 2 , in 
agreement with MC simulation. 

The decay is the only significant background source. It contributes either via in flight 
decays of the two pions {tt^ — > \i i^), or via a decay of a single pion and misidentification 
of another pion as a muon. Only pion decays resulting in muons almost collinear to the pion 
direction, and thus consistent with a three-track vertex and satisfying the total and transverse 
momentum requirements, contribute to background. 

Three methods of background evaluation are considered. 

• Owing to the symmetry properties of the detector, the kinematic distribution of the 
background events is to a good approximation identical to that of the reconstructed 
lepton number violating "wrong muon sign" (WS) 7r =F ^ ± // =l= candidates multiplied by 
a factor of 2. This observation allows to estimate the background contamination as 
(3.3 ± 0. 

5 s tat. i 0.5 sys t.)%, where the quoted statistical uncertainty is due to the limited 
number of data WS candidates (-/Vws = 52), and the systematic uncertainty is discussed 
below. 
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Figure 1: (a) Reconstructed spectrum of TT ± /j 1 + fi~ invariant mass: data (dots), MC sim- 

ulation and background estimate with WS events (filled areas), fit to background with an 
empirical function (smooth line), (b) Reconstructed spectrum of 3^ invariant mass: data (solid 
line histogram) and MC simulation (filled area). Signal regions are indicated with arrows. 



• MC simulation of the sample leads to a background estimate of (2.4 ± 0.7)%. This 
method also gives an estimate of the expected number of WS data events in the signal 
region: iV^g = 52.6 ± 19.8. The quoted uncertainties are systematic due to the limited 
precision of MC description of the high-mass region. They have been estimated from the 
level of data/MC agreement in the control reconstructed mass region of (465; 485) MeV/c 2 . 

• Fitting the mass spectrum in the region between 460 and 520 MeV/c 2 , excluding the signal 
region between 485 and 502 MeV/c 2 , with an empirical function (similar to that used in 
the E865 analysis [U [9]: a constant plus an exponentiated cubic polynomial) using the 
maximum likelihood estimator and assuming a Poisson probability density in each bin 
leads to a background estimate of 3.1%. 

The first method is considered the most reliable, and is used in the subsequent analysis. The 
background estimated with the WS and extrapolation methods is shown in Fig. [Ik.. The 
degree of agreement of the three background estimation methods is demonstrated in Fig. [2 

The distributions of the data and WS candidates in the {p\\ M^^) plane are dis- 

played in Fig. El The -p\ distributions of the data K wfJifJi candidates, MC events and 

background evaluated with WS events are shown in Fig. [U 

Normalisation sample 

The normalisation channel is well understood in terms of simulation, being of primary 
physics interest to NA48/2 |1 H I14| . The reconstructed 3^ invariant mass spectrum is presented 
in Fig. [lb: the non-Gaussian mass tails are due to ir^ — > i^v^ decays in flightH. The small 
deficit of MC events at low outside the signal region, explained in part by radiative 
decays, is not relevant for the present analysis. The background contamination is negligible. 
The number of candidates in the signal region is 2.386 x 10 9 which, taking into account 

1 Given that BR(if3 W )/BR(if^ MM ) ~ 10 6 , the filtering of the data stream and the analysis are performed is 
such a way that the candidates are effectively pre-scaled by a factor of 100, which strongly reduces the data 
volume. The shown mass plot is based on the pre-scaled sample. 
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Figure 2: Reconstructed M n ^ spectra of (a) K n ^ and (b) WS candidates: data (dots), 
and K % aa MC simulations (filled areas); fit to background using the empirical parameterization 
as explained in the text (solid line). The standard signal region is indicated with arrows. 



acceptance, trigger efficiency and BR |15j . corresponds to a number of kaon decays in the fiducial 
volume of Nk = 1-9 x 10 11 . The measured resolution of 03^ = 1.7 MeV/c 2 is in agreement 
with simulation, and is significantly smaller than a n ^ due to the smaller Q-value of the 
decay. 




3 Interpretation of the data 

Form factor parameterizations 

The decay is described as proceeding via single virtual photon exchange, resulting in a spectrum 
of the z = (M^^/Mk) 2 kinematic variable sensitive to the form factor W(z) [T]: 

where = tti^/Mk, r w = m n /MK, and A(a, b, c) = a 2 + b 2 + c 2 — 2ab — 2ac — 26c. The two- 
dimensional decay probability, which is used for evaluation of the geometrical acceptance with 
MC simulation, can be found for instance in Ref. [I]. The decay density is corrected by the 
Coulomb factor 

_ . , r 2naQ i Q j ( ^± ^ " J 

ij=l,2,3; i<j P ^ ^ J 

where Qi = ±1 are the electric charges of the daughter particles, < fyj < 1 are their relative 
velocities, and a is the fine structure constant. The relative velocities depend on invariant 
masses only; in particular, for the muon pair 0^ = 1 - [2r 2 /(z - 2r 2 )] 2 . The following models 
of the form factor W(.z) are considered. 

1. Linear: W(z) = GpM^foil + 5z) with free normalisation and slope (|/o|,<5). Decay rate 
and spectrum are not sensitive to the choice of the sign of /q. 



8 




Figure 3: Distributions of data events in the {pj,, M wfJjfJi ) plane for (a) K n ^^ candidates; (b) 
WS events. The signal region is enclosed within rectangles. The colour coding is in logarithmic 
scale. 
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Figure 4: Reconstructed p\ spectrum of candidates: data (dots), K wfJifM MC simulation 

and K^ w background estimated from WS events (filled areas). The upper limit of the signal 
region is indicated with an arrow. 
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Figure 5: (a) z spectrum of the selected candidates (dots) and K 3lT background estimated 

from WS events (filled area), (b) The reconstructed dT n ^/dz spectrum fitted to a linear form 
factor. 



2. Next-to-leading order ChPT [2]: W(z) = G F M^(a + + b + z) + W n7T (z) with free parameters 
(a + ,6 + ) and an explicitly calculated pion loop term W 7rw (z) given in [2]. 

3. Combined framework of ChPT and large- N c QCD (Ref. [3] and Appendix of Ref. |16|): 
the form factor is parameterized as W(z) = W(w, (3,z) with free parameters (w,/3). 

4. ChPT parameterization [4] involving meson form factors: W(z) = W(M a , M p , z). The 
resonance masses (M a , M p ) are treated as free parameters. 

Form factor measurement 

The reconstructed z spectrum of the data events is shown in Fig. [5^. The values of dT^^^/dz in 
each z bin, which can be directly compared to the theoretical expectations (pQ), are computed as 



Nj-N? _ A 3n (l-e 3w ) J_ ± 
N 3n Ai(l-Si) Azi t k 



(dT n ^/dz)i = l — * • ™\ ■ — ■ — ■ BR(K 3n ). (3) 



Here iVj and Nf are the numbers of K wflfl candidates and background events in the i-th bin, N 37T 
is the number of selected K 37T events, Ai (ranging from 5.2% to 25.2%) and £j are the geometrical 
acceptance and trigger inefficiency in the i-th bin for the signal sample, and ^t,- = 22.2% and £3,,- 
are those for K 3n events. Trigger inefficiencies are e ~ 0.2%, dominated by timing misalignment 
and therefore independent of kinematics and mostly cancelling in ([3]). The bin widths Azj for 
the {dT^^/ dz)i computation are chosen to be 0.02 for all bins except the first and last ones 
limited by the phase space: 4r^ < z < (1 — r n ) 2 . The resolution in the z variable is smaller than 
the bin width, ranging from at z = 4r^ to 0.007 at z = (1 — rv) 2 . The external inputs are the 
kaon lifetime tk and the branching ratio of the normalisation decay mode BR(A'3 7r ) [15J. 

The effective Z{ values, at which (dT^^^/dz)^ is evaluated, are corrected for non-linearity of 
dTnun/dz following |17j . This results in a significant correction in the first z bin where dT^^/dz 
has the largest gradient. 

The values of dT^nu/dz vs z and the result of the fit to the linear form factor are presented 
in Fig. [5b. The fits to the other models are very similar, and are not shown. 
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Figure 6: (a) Reconstructed cosOk^ spectrum of candidates: data (dots), MC 

simulation and background estimated from WS events (filled areas), (b) Data/MC ratio 
(background subtracted). The simulation does not include any forward-backward asymmetry. 



Model-independent branching fraction and asymmetries 

The mo del- independent BR is evaluated by integration of the spectrum ([3]) normalised to the 
full decay width H/tk- Separate measurements of the BR for K + and K~ decays allow 
the evaluation of the CP violating charge asymmetry of the decay rates: A(K^^) = (BR + — 
BR-)/(BR+ + BR-). 

Another interesting observable is the forward-backward asymmetry in terms of the angle 6k^ 
between the kaon and opposite sign muon three-momenta in the dimuon rest frame: 

A = iV(cos Bkjx > 0) ~ A(cos Oxy < 0) (4) 
FB N(caB0 Klt >0) + N(coB9 K p<Qy 

Non-zero values of the asymmetry can be obtained with the presence of scalar type interactions. 
Afb is expected to be sensitive to non-SM effects, and to be enhanced in the K nfJifJi decay with 
respect to the K^ ee decay [15] . 

Due to the limited size of the data sample, only the asymmetry integrated over the z vari- 
able is considered in the present analysis. The K + and K~ samples are summed up for this 
measurement. The background is subtracted using the WS events. The cosOk^ spectra 
(i.e. the sums of cosO K + ll - and cos0^- M + distributions) of the data, simulated signal 
and background events, and the data/MC ratio are displayed in Fig. [6] (no forward-backward 
asymmetry is simulated in MC). 

Systematic uncertainties 

The statistical errors include those due to the numbers of candidates and WS data events 
used to estimate the background. 

The scale factor of 2 applied to WS candidates to model the background in the sample 
has an uncertainty due to the dependence of Dalitz plot density and acceptance on kinematic 
variables. MC simulation of the decay mode predicts the scale factor to be (1.7 ± 0.3 s t a t.)- 
The stability of the ratio of spectra of data to WS candidates over the control 

region of (400; 480) MeV/c 2 suggests a similar uncertainty of 0.3. A dedicated study of the 
background to the — > ■K + ir~e ± v decay leads to a similar conclusion [19]. The effect of this 
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Table 1: Summary of the statistical and systematic uncertainties. 



Model 


Parameter 


Statistical 


Background 


Muon ID 


Pion ID 


External 


(1) 


l/ol 


0.039 


0.006 








0.002 




5 


0.56 


0.11 


0.01 








(2) 


a + 


0.038 


0.006 





0.002 


0.002 




b+ 


0.142 


0.027 


0.005 


0.006 


0.005 


(3) 


w 


0.014 


0.003 





0.001 









0.61 


0.12 


0.02 


0.02 


0.02 


(4) 


MjGeV 


0.083 


0.016 


0.002 


0.001 


0.001 




M b /GeV 


0.027 


0.005 


0.001 


0.001 


0.001 




BR x 10 s 


0.21 


0.07 


0.04 


0.08 


0.07 



uncertainty is quoted as the systematic error due to background subtraction. The background 
contamination is thus (3.3 ± 0.5 sta t. ± 0.5 sys t.)% = (3.3 ± 0.7)%. 

Muon identification efficiency measured from a — > \i^v sample collected during a special 
low intensity run with a loose 1-track trigger condition varies from 98.3% to 99.5%, depending on 
track momentum, in the signal momentum range above 10 GeV/c. It is reasonably well described 
by MC simulation: the difference between data and simulated efficiencies is below 0.4% in the 
whole analysis momentum range. This difference is propagated as the corresponding systematic 
uncertainty. 

Uncertainties due to the pion identification efficiency, which partially cancel due to sym- 
metrisation of the and selections, are evaluated by comparing the results with the 

pion identification criterion switched on and off for the K 3lT mode. 

The residual effects of the ~ 0.2% trigger inefficiency are negligible, as explained above. 
Uncertainties due to imperfect simulation of beam momentum and angular distributions are 
negligible, as these effects are small and largely cancel due to similarity of signal and normalisa- 
tion topologies. Sizable uncertainties arise from the external input: BR(i^3 7r ) = (5.59 ± 0.04)% 
is experimentally known with a limited relative precision of 0.7% |15j . 

Uncertainties on the model parameters and the model-independent BR are summarised in 
Table 1. Systematic uncertainties are significantly smaller than the statistical ones, owing to 
the limited size of the data sample. Systematic errors on the asymmetries A(K^ „) and Afb, 
which mostly cancel between the pairs of K^un, samples entering the asymmetry definitions, are 
negligible. 

Search for the lepton number violating decay K^ 1 — > Tv^fi^fi 1 * 1 

The decay — > 7r =F // =l= /i = ' = violating lepton number by two units can proceed via a neutrino 
exchange if the neutrino is a Majorana particle, and the E865 upper limit on its rate [9] currently 
provides the strongest constraint on the effective Majorana neutrino mass {m^) |201 12 1 J . This 
decay has also been studied in the context of super symmetric models with i?-parity violation |21j . 

A new upper limit on BR(iC ± — > 7r T /i ± /i ± ) can be established by analyzing the WS data 
mass spectrum shown in Fig. ^p. The expected background is estimated by MC simulation of the 
sample. The Feldman-Cousins method |22] is employed for confidence interval evaluation; 
the systematic uncertainty of the background estimate is taken into account. The geometrical 
acceptance is conservatively assumed to be the smallest of those averaged over the and 
K 3n samples (A n ^ = 15.4% and A 3n = 22.2%). 
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Table 2: The measured model parameters, their correlation coefficients, x 2 / n df of the fits, and 
the model-independent BR. 



Model (1) 




p(l/oU) = - 


-0.993 


X 2 /ndf = 12.0/15 








l/o| = 


0.470 


± 0.039 stat . 


± 0.006 syst . 


± 0.002 cxt . = 


0.470 


± 


0.040 


5 = 


3.11 


± 0.56 sta t. 


± o.n syst . 


= 


3.11 


± 


0.57 


Model (2) 




p(a+,b+) = 


-0.976 


X 2 /ndf = 14.8/15 








a+ = 


-0.575 


± 0.038 stat . 


± 0.006 syst . 


± 0.002 ext . = - 


0.575 


± 


0.039 


b+ = 


-0.813 


± 0.142 stat . 


± 0.028 syst . 


± 0.005 ext . = - 


0.813 


± 


0.145 


Model (3) 




p(w,P) = 0.999 


X 2 /ndf = 13.7/15 








w = 


0.064 


± o.oiw 


± 0.003 syst . 




0.064 


± 


0.014 


P = 


3.77 


± 0.6W 


± 0.12 syst 


± 0.02 ext . = 


3.77 


± 


0.62 


Model (4) 




P {M a ,M p ) -- 


= 0.999 


X 2 /ndf = 15.4/15 








M a /(GeV/c 2 ) = 


0.993 


± 0.083 stat . 


± 0.016 syst . 


± O.OOlext. = 


0.993 


± 


0.085 


M p /(GeV/c 2 ) = 


0.721 


± 0.027 stat . 


± 0.005 syst . 


± O.OOlext. = 


0.721 


± 


0.028 


BR x 10 s = 


9.62 


± 0.2W 


± O.llsyst. 


± 0.07 ext . - 


9.62 


± 


0.25 



Table 3: Summary of measurements of the form factor parameters with their statistical and 
systematic uncertainties. Some publications quote statistical errors only, see [16] for further 
details. 

Decay ^nee ^nee h ■ , up ^TTfifi 

Reference [7J P [3] p] [6] this result 



l/ol 


0.533 ±0.012 


0.531 ±0.016 


0.470 ± 0.040 


6 1.31 ± 0.48 


2.14 ±0.20 


2.32 ±0.18 2.45^:95 


3.11 ±0.57 


a + 


-0.587 ±0.010 


-0.578 ±0.016 


-0.575 ± 0.039 


b+ 


-0.655 ± 0.044 


-0.779 ± 0.066 


-0.813 ±0.145 


w 


0.045 ± 0.003 


0.057 ± 0.007 


0.064 ±0.014 


P 


2.8 ±0.1 


3.45 ± 0.30 


3.77 ±0.62 


M a /GeY 




0.974 ± 0.035 


0.993 ± 0.085 


Mp/GeV 




0.716 ±0.014 


0.721 ±0.028 



4 Results and discussion 

The measured values of the model parameters and the model-independent BR are presented in 
Table 2. The overall precision is limited mainly by the statistical uncertainties. Each of the form 
factor models provides a reasonable fit to the data, but the statistical precision is insufficient to 
distinguish between the models. The 68% confidence level contours for the pairs of parameters 
are presented in Fig. [71 overlayed with those obtained from the analysis of the NA48/2 K vee 
sample [16] . 

The measurements of the form factor parameters are summarized in Table 3. The present 
measurement is in agreement with the previous ones based on both [6] and K nee [7j [U [16] 

samples. Our measurement of the form factor linear slope 5 further confirms the contradiction 
of the data with the meson dominance models [23] which predict the slope parameter to be in 
the range from 0.5 to 0.9. The measurements of the a+ parameter (this result and the i^ee 
measurements [3 [16]) are in agreement with a theoretical expectation of a+ = — 0.6^ 'g [241 125] . 
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Figure 7: 68% confidence level contours for the parameters of the form factor models (larger 
contours) overlaid with those from the NA48/2 K nee analysis |16] (smaller contours). The 
dominant statistical errors only are considered. 

The measurements of BR(i^T 7r/J/J ) are summarized in Table 4. Our measurement is in agree- 
ment with those reported in [6j [10], but disagrees with the earliest one [5] by 4.5 standard 
deviations. It should be noted that the BR measurement [5], obtained under the assumption 
of the form factor slope measured with K nee [Z], is inconsistent with theoretical expectations 
based on K nee form factor measurements assuming that the K n u decays are dominated by the 
K — > 7T7* form factor [2]. The precision of our BR measurement represents a factor of ~ 3 
improvement with respect to the most precise earlier measurement [6j. 

The model-independent branching ratios measured separately for K + and K~ decays are 

BR+ = (9.70 ± 0.26) x 10~ 8 , BR~ = (9.49 ± 0.35) x 10~ 8 , 

where the quoted uncertainties are statistical only. Neglecting the systematic uncertainties of 
BR + and BR~ (which are small compared to the statistical uncertainties, and mostly common 
to the K + and K~ measurements), and the possible charge asymmetry of decay rates which 
is experimentally compatible with zero within 2 x 10~ 3 precision [26], we measure the charge 
asymmetry to be A(K^^) = (1.1 ± 2.3) x 10 -2 . This is a factor of ~ 5 improvement in precision 
with respect to the only previous measurement [10] . and is compatible with CP conservation. A 
limit for the charge asymmetry of |A(ir^ )| < 2.9 x 10~ 2 at 90% CL can be deduced from the 
above value. The experimental precision is far from the SM expectation \A(K± lfl )\ ~ 10" 4 [2] 
and even the SUSY upper limit \A(K± )\ ~ 10" 3 [23 [28] for the CP violating asymmetry. 
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Table 4: Summary of BR(i^ 7r/i/1 ) 


measurements. 


Reference 


Beams 


K^nn candidates 


BR x 10 8 


E787 [SJ 


K+ 


207 


5.0 ±0.4 ±0.6 ±0.7 


E865 [6] 


K+ 


430 


9.22 ±0.60 ±0.49 


HyperCP [TO] 


K± 


110 


9.8 ± 1.0 ±0.5 


NA48/2 


K± 


3120 


9.62 ±0.21 ±0.11 ±0.07 



The forward-backward asymmetry has been measured to be Apg = (—2.4 ± 1.8) x 10 2 , 
where the error is dominated by the statistical uncertainty. It corresponds to an upper limit 
of \A FB \ < 2.3 x 10~ 2 at 90% CL. The achieved precision does not reach the upper limits for 
the SM contribution via the two-photon intermediate state if ± — > 7r =l= 7*7* — > tt^^ n~ |29j and 
MSSM contribution [18], which are both of the order of 10 -3 . 

Finally, -/Vws = 52 WS data events are observed in the signal region with the expected 
background of = (52.6 ± 19.8) estimated from MC simulation. Conservatively assuming 

the expected background to be 52.6 — 19.8 = 32.8 events to take into account its uncertainty, 
this translates into an upper limit of 32.2 signal events at 90% CL, leading to an upper limit of 
BR(if ^ — > tt^ n^ji^ ) < 1.1 x 10~ 9 at 90% CL. This is an improvement by almost a factor of 3 
with respect to the best previous limit [9], allowing a bound on the effective Majorana neutrino 
mass of (m w ) < 300 GeV/c 2 to be established [20] . 

Conclusions 

From a sample of 3120 — > tt ± /j 1 + /j,^ decay candidates with (3.3 ± 0.7)% background con- 
tamination, the model-independent branching fraction has been measured to be BR = (9.62 ± 
0.21 s tat. ±0.11 sy st. ±0.07 ex t.) x 10~ 8 , and the form factor which characterizes the decay has been 
evaluated in the framework of four models. Upper limits for the CP violating charge asym- 
metry and (for the first time) the forward-backward asymmetry of the decay rate have been 
established. An upper limit of 1.1 x 10~ 9 for the branching fraction of the lepton number violat- 
ing — > 7r =F ^ ± // ± decay has been obtained. The achieved precisions dominate the currently 
available measurements. 
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